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ABSTRACT 

Gravitational, magnetic and superfluid forces can stress the crust of an 
evolving neutron star. Fracture of the crust under these stresses could affect 
the star's spin evolution and generate high-energy emission. We study the 
growth of strain in the crust of a spinning down, magnetized neutron star and 
examine the initiation of crust cracking (a starquake) . In preliminary work 
(Link, Franco & Epstein 1998), we studied a homogeneous model of a neutron 
star. Here we extend this work by considering a more realistic model of a 
solid, homogeneous crust afloat on a liquid core. In the limits of astrophysical 
interest, our new results qualitatively agree with those from the simpler model: 
the stellar crust fractures under shear stress at the rotational equator, matter 
moves to higher latitudes and the star's oblateness is reduced. Magnetic stresses 
favor faults directed toward the magnetic poles. Thus our previous conclusions 
concerning the star's spin response still hold; namely, asymmetric redistribution 
of matter excites damped precession which could ultimately lead to an increase 
in the spin-down torque. Starquakes associated with glitches could explain the 
permanent offsets in period derivative observed to follow glitches in at least 
three pulsars. 

Subject headings: dense matter — magnetic fields — stars: magnetic fields - 
stars: neutron — pulsars: individual (Crab) 
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1. Introduction 

Stresses developing in the crust of an evolving neutron star could fracture the crust 
(i.e. a starquake), possibly affecting the star's spin evolution and generating high-energy 
emission. Spin down (e.g., in isolated pulsars; Paym fc Pines 1971| ) or spin up (e.g., in 
accreting neutron stars) changes the equilibrium shape of the star, building stresses. In 
"magnetars", decay of the superstrong field (B> 10 14 G) could break the crust and produce 
episodes of intense gamma-ray emission ( [l'hompson fc Duncan 1996| ; [1'hompson fc Blacs| 



|1998[ ). Differential rotation between the crust and the interior neutron superfluid may also 
stress the crust to its breaking point ([Ruderman 1976| ). 



Starquake-induced rearrangement of the star's mass distribution changes the star's 
moment of inertia causing precession and polar wandering (Link, Franco & Epstein 1998; 
hereafter LFE). Anomalous spin behavior might therefore signal the occurrence of a 
starquake. Such evidence for crust cracking may already exist in several isolated pulsars. 
Permanent offsets in period derivative following glitches have been observed in the Crab 
pulsar flLyne, Pritchard, fc Smith 1993|) , PSR1830-08 and probably PSR0355+54 ( (Shemar 



|fc Lyne 1996|) . The observed offsets all have the same sign, and correspond to increases in 



spin-down rate. In LFE we interpreted these offsets as due to permanent increases in the 
torque acting on the neutron star. We showed that cracking and readjustment of the stellar 
crust in response to the star's slow down can increase the angle between the star's spin and 
magnetic axes, leading to an increase in the spin-down torque in some models of pulsar spin 
down. 

Other evidence for cracking of the neutron star crust can be found in the Soft Gamma 
Repeaters, thought to be strongly magnetized neutron stars (.B^IO 14 G), or magnetars. 
Unlike radio pulsars which have weaker magnetic fields and are powered by rotation, 
magnetars are thought to be powered by their intense magnetic fields. Duncan & Thompson 
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(1994; also Thompson & Duncan 1995; 1996) have suggested that SGR outbursts represent 
cracking of the crust by magnetic stresses. Evidence in favor of this hypothesis is the 
striking statistical similarities between bursts in SGRs and earthquakes QCheng et al. 



1996| ; Gogus et al. 1999| ). For both phenomena the energy released per event obeys similar 
power-law scalings and the waiting times between events are strongly correlated with one 
another. 

The potential importance of violent crust dynamics in neutron star evolution motivates 
this detailed study of spin-down induced starquakes. The equilibrium shape for a spinning 
star is an oblate spheroid which becomes more spherical as the star slows down. The liquid 
core is able to change its shape smoothly as the star slows down. In contrast, the solid 
crust is strained as it evolves through a continuous series of equilbrium configurations, and 



eventually will crack if it is brittle. In [LFEj , we modeled the neutron star as a homogeneous, 
self-gravitating, elastic sphere and followed the evolution of strain in its crust as the star 
spins down. From our calculations we developed the following picture of starquakes. Crust 
cracking occurs as equatorial material shears under the compressive forces arising from the 
star's decreasing circumference. The star's oblateness suddenly decreases as matter moves 
to higher latitudes along a slip or fault plane that is perpendicular to the stellar surface and 
crosses the equator at an angle of ~ 30 — 45° (see Fig. 6). Since magnetic stresses suppress 
shearing near the magnetic poles and across the field lines, starquakes likely originate near 
the two points on the equator farthest from the magnetic poles and propagate preferentially 
toward the magnetic poles. The matter redistribution breaks the rotational symmetry 
of the star, causing the star to precess. Damping of the precession eventually restores 
alignment between the angular velocity and angular momentum, ultimately increasing the 
angle between the spin axis and the magnetic moment. In some models of pulsar spin-down 
{e.g., the vacuum-dipole model) an increase in the alignment angle between the rotation 
and magnetic axes increases the spin-down torque. The magnitude of the starquake-induced 
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changes in the alignment angle are consistent with the changes in the spin-down rates seen 
in the Crab pulsar ( |LFE|) . The observations of spin-down offsets in coincidence with glitches 
indicates some connection between structural readjustments and glitches. For a recent 
discussion of this possibility, see Link & Epstein (1996). 

In this paper we extend the analysis to the more realistic case of a brittle shell floating 
on a liquid core. We find that the stresses induced by the star's spin down produce similar 
behavior to that found in LFE. Crust cracking associated with other sources of stress, such 
as magnetic or superfluid stresses, will be considered in future work. The organization of 
this paper is as follows. Section |2| presents a qualitative description of the growth and 
release of strain in the crust of an idealized, spinning-down neutron star. In we calculate 
the response of the crust's shape as it spins down, treating the star as a homogeneous solid 
afloat on a liquid core. In §[|, we explore the cracking and matter flow for a model with 
realistic crust thickness. In §[| we discuss the effects of a magnetic field anchored in the 
crust. Section || summarizes our results. 



2. Growth and Release of Strain in a Spinning-down Star 

Our analysis is based on the assumption that the neutron star crust is brittle, and 
therefore cracks when sufficiently strained^. We consider a neutron star spinning with 



2 Experience with terrestrial materials shows that material properties at high pressure 
can differ from those at low pressure. For example, terrestrial materials typically become 
ductile when subjected to pressures in excess of their shear moduli. Nevertheless, deep-focus 
earthquakes are known to originate from regions of very high pressure flGreen fc Houston 



1995|) . These faults are thought to be triggered by densification phase changes; small regions 



of higher density nucleate as the material is stressed, and act as a lubricant for shearing 
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angular velocity Q which undergoes a change — SQ in its rotation rate where < 5Q <C Q. 
Our goal is to describe the effects on the star's crust resulting from this change in rotation 
rate. In equilibrium, the star has a spheroidal shape with an equatorial bulge of relative size 
~ (QR/vk) 2 where R is the radius of the non-rotating configuration and vu is the Keplerian 
velocity. As the star spins down, the fluid interior readjusts its shape continuously, following 
(in the homogeneous case) a sequence of Maclaurin spheroids (see, e.g., Chandrasekhar, 
S. 1987; p. 77). As the solid crust attempts to readjust its shape, however, strain builds. 
Material originally at r is displaced to r + u(r), where u(r) is the displacement field. The 
local distortion of the solid is given by the strain tensor (see, e.g., Landau & Lifshitz 1959; 

In a local coordinate system in which this matrix is diagonal, the eigenvalues en represent 
compression (en < 0) or dilation (en > 0) along the respective axes. Let the largest and 
smallest eigenvalues be e\ and e s , respectively. The strain angle is e\ — e s , and we refer to 
the plane containing the corresponding principal axes as the stress plane. In Fig. 1 we show 
schematically a block of matter that is compressed along the y direction and allowed to 
expand along the z direction. The stress plane then is parallel to the y — z plane. 

The elastic limit of the matter is reached first in the regions of the crust where the 
strain angle is maximum. There the material cracks, forming a fault plane perpendicular 
to the stress plane (Fig. 1). Matter slips along the fault plane and redistributes in a way 
that reduces the equatorial circumference, thus releasing the strain. The fracture geometry 
depends on where the strain angle is maximum and how its value changes over the crust. 

motion. Analogous processes might occur in the high-pressure material of the neutron star. 
E. Ramirez- Ruiz and R. I. Epstein are examining whether the phase transition from spherical 
nuclei to rod-like nuclei may facilitate starquakes. 
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We now turn to a calculation of the spin-down induced strain in the neutron star crust. 

3. Matter Redistribution in a Spinning-down Star 

To explore the accumulation and release of strain in a slowing neutron star, we model 
the star as a two-component homogeneous spheroid; a self-gravitating core of incompressible 
liquid on which rests a crust of uniform density. The non-rotating configuration is a sphere 
of core radius R' and total radius R. 

3.1. Equilibrium Equations 

A spinning, self-gravitating object has an equilibrium configuration given by 

V • T - pV0 = F, (2) 

where T is the material stress tensor, <ft is the gravitational potential per unit mass and 
F is the centrifugal force density associated with the rotation. For rigid material, strain 
develops as the force is changed. The strain in the material is given by the strain tensor of 
eq. ([]]). For incompressible matter we have 

V • u = uu = 0, (3) 

where repeated indices indicate sums. 

For sufficiently small strains, the material can be described using Hooke's Law in terms 
of the stress tensor (see, e.g., Landau & Lifshitz 1959, p. 11). 

Jij = -p5ij + an (4) 
Oij = 2/iUij, (5) 
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where p is the isotropic pressure, \i is the elastic shear modulus, represents the 
contribution of the material's rigidity to the stress, and we have applied eq. (|3|). 

We consider an initial configuration in which the material is under great internal 
pressure but unstrained (uij = 0). The equation of hydrostatic equilibrium is then 

- Vp - pV0 o = F , (6) 

and, in the limit of slow rotation, RQ/vk <C 1, F is small compared to both Vpo and pV0o- 

To derive the condition for hydrostatic equilibrium in the new configuration produced 
by the change in F, we consider a Lagrangian perturbation of the initial state. A 
Lagrangian perturbation A is evaluated moving with the matter, and is related to an 
Eulerian perturbation S through 

A = 5 + u-V. (7) 

A Lagrangian perturbation commutes with spatial derivatives as (see, e.g., Shapiro & 
Teukolsky 1983; p. 131) 

A ® — ® A ^ Uj V (8) 
dxi dxi dxi 3 

We obtain the new equilibrium by taking a Lagrangian perturbation of eq. (H), i.e., 

A [V • T — pV0 = F] . (9) 
Using eqs. ([p) through @, we obtain, 

- VAp - pc 2 t V x V x u - P W5<P = 5F- V(u ■ Vp ) + (u ■ V)F - u k — (10) 

OXi 



where we have neglected terms of order u 2 and c t = yp/p is the transverse sound 
speed, which we take to be constant throughout the solid. For the case of slow rotation, 
RQ/vk 1, the last two terms in eq. ( P!U[ ) are smaller than the leading terms by a factor 
~ (QR/vk) 2 , and we neglect them. 
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We express the Eulerian change in centrifugal force by means of a scalar potential \ 
defined as 

6F = pV X . (11) 

Introducing a second potential h, 

c]h = --Ap-50-x+ -u- Vp , (12) 
P P 

we rewrite eq. fllPP as 

V x V x u = Vh. (13) 
Taking the divergence of this equation, we see that the potential satisfies 

V 2 /i = 0. (14) 



The displacement u is obtained by solving eqs. (0) and (|14D with ([121) subject to boundary 
conditions. 



3.2. Displacement Field of the Crust 

In a spherical coordinate system (r, 9) with Q directed along the z-axis, the centrifugal 
potential associated with a change in rotation rate 5Q is 

X = l^Qr 2 (l-P 2 (9)), (15) 

where P 2 = (3cos 2 # — l)/2 is the second Legendre polynomial. The P 2 (0) symmetry of \ 
suggests a solution to eq. ( Ill ) of the form 



h(r,9) = ^Ar 2 + ^P 2 (6) + C, (16) 

where A, B and C are constants. For the displacement field, we seek a solution of the form 

u r (r,6) = f(r)P 2 (6) (17) 
u 6 (r,6)=g(r)^P- (18) 
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where /(r) and g(r) are functions of depth in the star only. The solution to eqs. ( p"3"D and 
(0), with (|16D, and subject to condition (|3|), is 

u r (r, 6) = (ar - yr 3 - JL + A) p a(e) (ig) 

where a and 6 are additional constants. We require four boundary conditions to specify the 
four coefficients. 



3.3. Boundary Conditions 

At the outer boundary the crust encounters a vacuum, while at the inner boundary 
it is in contact with the liquid core. Since neither of these media support traction, two 
boundary conditions follow from requiring that no shear exist at each boundary, i.e., 

cx r9 = at r = R, R', (21) 



or, 



due Ua 1 du r . , , . 

- — + = at r — R, R. (22) 



dr r r 89 
In terms of the coefficients 

8 t , B 8 b 

a AR 2 + = 23 

21 2R 3 3 R 5 K J 

8 . „,o B 8 b , 

a AR' 2 + - = 0. 24 

21 2R' 3 3R' 5 K J 

Two other boundary conditions follow from requiring that the material be in mechanical 
equilibrium at the boundaries. This requirement is satisfied by ensuring that 
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is continuous at r = R, R'. Hence, 



n 9 du r 



J r=R-e 



„ 9 9« r 



[-Ap] r= ^ 



r=R'-e ' 



(26) 



(27) 



J r=_R'+e 

where e is an infinitesimal positive distance. We have required Ap = just above the star 
and have taken c t = in the liquid core. These jump conditions specify the material stress 
required to maintain a pressure difference across a boundary. 

The local pressure change follows from eq. ( |T2"D 



Ap = —pc t h — p5(p — px + u • Vpo- 



(28) 



The initial, unstrained configuration is that of a Maclaurin spheroid whose shape, for small 
eccentricity e, can be written 

1 



R{6) = R 



l--e 2 P 2 (6) 



= R + e(6). 



(29) 



The gravitational potential for a slowly-rotating Maclaurin spheroid is (Shapiro & Teukolsky 
1983, Black Holes, White Dwarfs and Neutron Stars, p. 169) 



2R 2 



-r z + 



4r s 



■e(6>) 



(30) 



3 5R 

Since the crust of a neutron star contains only ~ 1% of the star's mass, we neglect the 
crust's effect on the gravitational potential. As the star spins down, e changes by u r (R,6). 
Hence, the Eulerian change in the potential is 



(31) 



4 r 2 

"=^Gp—u r (R,8). 

5 R 



For slow rotation, Vpo is nearly radial, and we may approximate the last term in eq. (p8[) as 

4 

u ■ Vp ~ -pu r {r, 9)V(po ^ --irGp 2 ru r (r, 6), 



(32) 
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where we used the spherical (non-rotating) gravitational potential in the last step, which is 
justified in the limit of slow rotation. 



Combining eqs. 
the surface yields 



, (|32|), and (|T^) through (|TTD, the boundary condition (^) at 



2v 2 k f(R) 2^M1R 



+ 



5 c 2 R 3 c 2 



- 2f(R) - 

O Ct 

where f(R) is defined by eq. ([H|) and v\ = AtcGR 2 p/3. 



AR + 



B 



(33) 



To match the solutions at the crust-core boundary, we require an expression for the 
pressure change in the core. We obtain it from eq. (|T0|) setting c t = since a liquid cannot 
support shear. We have 



5(f) + -u • Vpo - x 

P P 



0. 



Hence, 



Ap 
P 



H — u ■ Vpo — X = constant 
P 



(34) 



(35) 



in the core. The pressure change in the shell obeys 



Ap 



5(p + -u • Vp - x = (k 



A 2 B 

Ar 2 + — 



P 2 (0) + C. 



P P 

Thus, the boundary condition (|27| ) at the crust-core boundary, using eqs. (|35| ) and (|36l) 
gives 



(36) 



du r 
dr 



AR' 2 + 



B 

^3 



Pi{6) + constant. 



(37) 



r=R'+e 

^From the form assumed for u r , we see that the constant is zero. We finally obtain 



f'(R') = -\ 



(38) 



Eqs. (0), (P0|) , and the boundary conditions (0), fl24|) , (|33D and (pq) complete the 
description of the displacement field driven by a change in rotation rate — 5Q. 
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Recent calculations with realistic equations of state suggest that neutron stars have 
crusts of thickness 0.05 - 0.08i? ( [Akmal et al. 1998| ; |Lorenz et al. 1993] ); in the following, we 
take t = R — B! = 0.05i?. The corresponding displacement field is shown in Fig. 2. From 
the displacement field we can calculate the strain distribution resulting from spin down. 
The strain tensor ([!]) can be diagonalized to obtain the strain angle and the planes along 
which the material can shear. The geometry of crust cracking can be obtained from an 
analysis of the distribution of strain angles over the crust, to which we now turn. 



4. Cracking the Crust 

As the star spins down and the strain increases, the crust breaks when the local strain 
angle reaches a critical value. Fig. 3 shows contours of constant strain angle through the 
crust; the strain angle is maximum at the crust-core boundary, in the equatorial plane. Fig. 
4 shows the three eigenvalues of the strain tensor at the inner boundary of the crust plotted 
as a function of latitude. The strain angle, shown by the darker curve, shows discontinuous 
derivatives where the maximum and minimum eigenvalues change. The dashed curve 
corresponds to the strain angle for a homogeneous elastic sphere model ( |LFri| ) at the same 



depth for comparison. In the equatorial plane (8 = 7r/2), the stress plane is the #0-plane. 
An element of matter on the equatorial plane is compressed by adjacent elements on the 
equator and expands in the r and 9 directions. Material compressed in this way shears 
along one of two planes that takes an angle of 7 ~ 30 — 45° with respect to the compressive 
force QGreen fc Houston 1995| ; see Fig. 1). Once the crust breaks, matter moves along 



faults F or F', inclined at an angle 7 to the rotational equator (see Fig. 6). Since the strain 
field has azimuthal symmetry, the break is equally likely to begin at all points on a ring 
corresponding to the intersection of the crust core boundary with the equatorial plane. The 
inclusion of magnetic stresses, discussed below, breaks this symmetry. 
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The propagation of the crack is a nonlinear process; cracking affects the strain field, 
which in turn affects the propagation of the crack. In general, cracking increases the stress 



on unbroken neighboring material QHertzberg 1996Q , causing the fault plane to grow in 



both width and length. Hence, the crack begins at the base of the crust in the equatorial 
plane, propagates toward the surface, and continues to grow up to some latitude. When the 
cracking ends, the strain field is everywhere subcritical. The width and length of the fault 
plane cannot be determined from the linear analysis given here. 

In this treatment, we have ignored the fact that the neutron star crust is stratified. 
Inspection of eq. [ID] shows that the most important material property controlling the 



development of strain is q; in particular, where q is small, the strain must be large. 
Strohmayer (et al. 1991) calculated c t as a function of density in the neutron star crust 
and found that it varies by a factor of ~ 10 2 . Cracking of the crust is favored in regions 
where c± is small, possibly causing the quake epicenter to be somewhat higher than the base 
of the crust. We expect that the effects of density stratification will not affect the overall 
geometry of crust cracking described here. 



5. Magnetic stresses 

We now examine the effects of the magnetic field that is anchored to the crust of the 
neutron star. The displacements induced by the spin down distort the magnetic field, 
generating magnetic stresses. If the magnetic field is not precisely aligned with the rotation 
axis, the magnetic stresses will break the azimuthal symmetry of the strain distribution in 
the crust. 

A complete description of the strain field that develops in magnetized material is 
obtained by solving the elasticity equations and Maxwell's equations self-consistently, 
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subject to boundary conditions at the star's surface. The surface boundary conditions, 
however, are complex to handle; in a realistic neutron star they depend on the exact form 
of the density gradients near the surface. To qualitatively illustrate the influence of the 
star's magnetic field on where starquakes originate and how they propagate, we consider 
the simpler problem of the magnetic fields effect on the material strain in an unbounded 
medium. 

In the absence of the field, the equilibrium state is given by 



where Fj represents gravitational and centrifugal forces. The presence of a magnetic field 
introduces several forces. First, the unperturbed field can exert a force on the crust. For 
this study, we assume that the initial field B is dipolar and hence force free; that is, 
the unperturbed Maxwell stress tensor Tjj obeys dTij/dxi = in the crust. Next, the 
displacement field in the crust distorts the field by <5B, generating a perturbation STij to 
the Maxwell stress tensor .0 The new equilibrium is given by 



where ba^ is the magnetically induce change in the material stress tensor. For an infinite 
medium (or a finite one in which boundary effects are unimportant), the solution is 



Magnetic stresses modify the displacement u by Su. For a magnetic field that is 
sufficiently small that this correction to the displacement field is small (\5u\ <C |u|), the 

3 The ratio of magnetic to material stresses is ~ f3 = (tu/c*) 2 where va is the 
Alfven speed and q is the transverse sound speed. In the outer crust, we have j3 = 
8 x 10" 4 (B/10 12 G) 2 (p/10 10 g cm" 3 )' 1 (q/10 8 m s _1 )~ 2 and we are justified in treating the 
magnetic effects perturbatively. 




(40) 



8 an 
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perturbation to the Maxwell stress tensor is 

SJ ij = T- i SB i B j + B i SB j ~ BtSB^] , (41) 

with 

SB = V x (u x B). (42) 

Here u is the displacement field in the absence of the magnetic field. For incompressible 
matter, the correction Suij to the strain tensor is 

Suij = — (day - ^ a iij ' ( 43 ) 

where /i is the shear modulus. The eigenvalues of this tensor give the correction to the 
strain angle due to the induced magnetic forces. Using the displacement field of eqs. (IT 



and (^(j) in eq. (f42|) , we obtain an estimate of the corrections to the strain angle on the 



star's rotational equator; these are shown in Fig. 5. 

The magnetic field affects our picture of crust cracking in two ways. First, the existence 
of the magnetic field strengthens the material, making it harder to shear. The material is 
weaker and hence most likely to shear where the magnetic field is weakest. Consider a star 
with a dipolar magnetic field inclined at an angle a to the rotation axis. In this case the 
material may begin shearing at the two points on the rotational equator that are farthest 
away from the magnetic poles as shown in the schematic of Fig. 6. Second, as the material 
attempts to move along the fault planes, the magnetic field opposes displacements that 
cross field lines and hence favors matter motion long fault F over fault F 1 . This motion 
moves matter closer to the magnetic poles. This picture is essentially the same as found for 
a simpler homogeneous sphere model (|LFE|) . 
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6. Discussion 

We have examined the evolution of strain in the crust of an idealized spinning-down 
neutron star and the initiation of starquakes as the material reaches critical strain. We 
modeled the neutron star as a self-gravitating core of incompressible liquid that supports 
a brittle crust. Crust cracking occurs as material shears under the compressive forces 
arising from the star's decreasing circumference. The star initially cracks near the star's 
rotational equator, at the base of the crust. The shearing motion along the fault decreases 
the oblateness of the star and pushes matter to higher latitudes. Magnetic stresses suppress 
shearing near the magnetic poles and across the field lines. Starquakes thus originate near 
the two points on the equator farthest from the magnetic poles and propagate toward 
the magnetic poles. Our conclusion is qualitatively the same as found in our simpler 
homogeneous spherical model ( [LFE ) . 

As material slides along a fault directed toward the magnetic poles, the star's principal 
axis of inertia shifts away from the magnetic axis, exciting precession. Eventually, damping 
restores alignment between the principal axis of inertia and the angular momentum axis, 
increasing the angle between the rotation and magnetic axes. In some models of pulsar 
spin-down, e.g., the magnetic dipole model, this growth in the alignment angle increases 
the spin-down torque. We might, therefore, expect to see a tendency of alignment angle to 
increase with age in young pulsars. In older pulsars, however, Tauris & Manchester (1998), 
find no relationship between magnetic alignment angle and spin-down torque. 

One possible observational signature from a starquake is a change of the pulse profile. 
A jump in the alignment angle a changes the duration of the line of sight's traverse through 
the pulse emission cone. In the Crab pulsar, the relative change in pulsed flux could be of 
order 1% (Link & Epstein 1997). Additionally, the energy released in a starquake could 
produce an observable enhancement of the star's luminosity. The energy release is of 
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order E ~ 10 42 (/i/10 31 erg cm~ 3 )(<i/10 5 cm) 3 (0 c /lO -2 ) 2 ergs, where /j, is the average shear 
modulus along a fault of length d and 9 C is the strain angle at which fracture occurs. Some 
fraction of the seismic energy will be damped, heating the crust. Eventually, a soft x-ray 
thermal wave will emerge from the star's surface. For localized deposition of 10 42 ergs at a 
density of 10 14 g cur 3 , Tang (1999) calculates a 5% enhancement of the star's luminosity 
occurring about 3 yr after the quake for a star with an internal temperature of 10 7 K and a 
stiff equation of state. For deposition at 10 13 g cm" 3 in a star with an internal temperature 
of 10 8 K, the enhancement is ~ 22%. The Chandra X-ray Observatory should be capable 
of detecting this enhancement in a closer source, such as the Vela pulsar. 

A starquake could drive non-thermal emission as well. Seismic energy that reaches the 
stellar surface couples to Alfven waves which propagate into the magnetosphere. Conversion 
of Alfven waves into 7-rays could occur if the waves are charge-starved or if their amplitudes 
are comparable to the background field strength of the magnetosphere (Blaes et al. 1989). 

We thank A. Olinto for valuable discussions and C. Miller for a critical reading of the 
manuscript. This work was performed under the auspices of the U.S. Department of Energy, 
and was supported in part by NASA EPSCoR Grant #291748, NASA ATP grant # NAG 
53688, by IGPP at LANL and by the Center for Thermonuclear Flashes at the University 
of Chicago. 
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stress plane 



Fig. 1. — Breaking of a compressed matter element. A block of matter subjected to horizontal 
compression and vertical tension shears along a fault plane as shown when critical strain is 
reached. Shearing along a complementary plane flipped over with respect to the plane shown 
is equally likely for isotropic material. The y-z plane is the stress plane. 
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Fig. 2. — The displacement field in a spinning-down neutron star. A cross 
the center of the star is shown for a model with relative crust thickness of 5% 
diameter decreases while the polar diameter increases. 
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Fig. 3. — Contours of constant strain angle as a function of depth in the crust for a non- 
magnetized neutron star (arbitrary units). The maximum strain angle occurs at the equator 
and on the inner boundary of the stellar crust. 




Fig. 4. — Strain eigenvalues and strain angle vs. latitude at the stellar surface for B = 
(arbitrary units). The curves show the strain angle and the eigenvalues as marked. The 
dashed line corresponds to the strain angle for the homogeneous sphere of LFE. 
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Fig. 5. — The strain angle of surface equatorial material in the presence of a magnetic field 
(arbitrary units). The curves correspond to different values of the angle a between the 
magnetic and rotation axes. The magnetic poles are at azimuthal angles n/2 and 37r/2. 
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Fig. 6. — A starquake. In the absence of a magnetic field, the material is equally likely 
to begin breaking along faults F and F', anywhere on the equator. These fault planes are 
perpendicular to the stellar surface. In the presence of magnetic stresses, fault F is more 
likely, creating "mountains" (indicated by the snow-capped peaks) and shifting the largest 
principal axis of inertia to a new direction e\ (fixed in the star). 



